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ABSTRACT 


This  study  examines  enlisted  attrition  behavior  for  the 
U.S.  Coast  Guard  and  develops  a  model  that  projects  attrition 
figures.  Survival  analysis  techniques  are  used  to  analyze  the 
empirical  attrition  behaviors  associated  with  an  individuals 
sex,  race,  marital  status,  and  military  occupational  skill 
(MOS)  .  In  this  study  males  tend  to  have  higher  survival 
probabilities  than  females,  non- Caucasians  higher  than 
Caucasians,  and  married  persons  higher  than  those  not  married. 
Aviation  MOSs  have  the  highest  survival  probabilities  and 
technical  MOSs  have  the  lowest. 

Modelled  survivor  functions  are  developed  for  two 
paygrades  because  each  contain  small  personnel  Inventories. 
These  modelled  survivor  functions  do  not  fit  the  data  as  well 
as  desired  but  are  nonetheless  used  pending  the  development  of 
sharper  alternatives. 

Finally  a  counting  model  based  on  the  Binomial 
Distribution  is  developed  that  projects  monthly  enlisted 
attrition  figures. 
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I .  INTRODUCTION 


A.  PROBLEM  STATEMENT 

The  United  States  Coast  Guard  (CG)  Is  a  military  service 
operating  within  the  Department  of  Transportation.  Its  main 
missions  are  aids  to  navigation,  boating  safety,  defense 
operations,  maritime  law  enforcement ,  environmental  response, 
marine  inspection,  ice  operations,  marine  licensing,  marine 
science,  search  and  rescue,  port  security  and  safety,  and 
waterways  management.  The  CG  maintains  an  active  duty 
hierarchical  personnel  structure  with  approximately  6,000 
officers  and  30,000  enlisted  members.  The  CG  Personnel 
Workiforce  Planning  (PWP)  office  located  in  Washington,  D.C. 
has  the  overall  responsibility  for  projecting  the  active  duty 
enlisted  personnel  needs  of  the  CG.  Their  responsibilities 
include  the  calculating  of  personnel  stocks,  promotion 
requirements,  attritions,  and  recruitment  needs.  The  current 
method  for  predicting  monthly  enlisted  attrition  by  paygrade 
is  to  calculate  a  running  average  or  mean  of  past  monthly 
attrition  and  use  this  mean  as  a  forecast  of  future  attrition. 
This  forecasting  method  is  crude,  ignores  many  issues,  and  is 
not  accurate  enough  to  meet  the  needs  of  the  CG. 
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Why  does  a  person  separate  from  the  CG?  What  are  the 
characteristics  associated  with  this  person?  Are  these 
characteristics  the  saune  for  each  paygrade?  Previous  CG 
studies  have  focussed  on  the  first  question  but  have  not 
looked  at  attrition  data  analytically.  The  first  question 
will  not  be  addressed.  This  thesis  will  look  at  the  next  two 
questions  and  will  propose  a  method,  based  on  a  survival 
analysis  of  CG  enlisted  attrition  data,  for  calculating 
monthly  personnel  attrition.  This  study  is  a  first  attempt  at 
examining  CG  enlisted  attrition  using  data  analysis  tools.  It 
is  designed  to  provide  a  procedure  for  studying  attrition,  and 
for  developing  a  basic  model  that  calculates  monthly  attrition 
which  can  be  expanded  by  future  studies. 

B.  COAST  QCARD  ENLISTED  PERSONNEL  BACKGROUND  INFORMATION 
The  Coast  Guard  enlisted  personnel  structure  is  composed 
of  9  paygrades,  E-1  to  E-9,  lowest  to  highest  respectively. 
(There  is  an  E-10  position  filled  by  one  person  but  this 
thesis  does  not  address  this  paygrade.)  Paygrade  E-1  is 

composed  of  persons  attending  CG  basic  training  school  and 
paygrade  E-2  is  composed  of  persons  who  have  completed  basic 
training  school,  have  been  assigned  to  active  duty  commands, 
and  are  preparing  to  attend  a  CG  specialty  school  (also  called 
A- school) .  Paygrades  E-S  through  E-9  are  divided  into  24 
active  duty  military  occupational  skills  (MOS)  or  sub¬ 
specialties  listed  below: 
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•  Aviation  Machinist's  Mate  (AD) 

'  Aviation  Electrician's  Mate  (AE) 

-  Aviation  Structural  Mechanic  (AM) 

-  Aviation  Survlvalman  (ASM) 

-  Aviation  Electronic  Technician  (AT) 

-  Boatswain's  Mate  (BM) 

-  Deunage  Controlman  (DC) 

>  Electrician's  Mate  (EM) 

-  Electronics  Technician  (ET) 

-  Fire  Control  Technician  (FT) 

-  Gunner's  Mate  (GM) 

-  Health  Services  Technician  (HS) 

-  Machinery  Technician  (MK) 

-  Musician  (MU) 

-  Marine  Science  Technician  (MST) 

-  Public  Affairs  Specialist  (PA) 

-  Quartermaster  (QM) 

-  Radarman  (RD) 

-  Radioman  (RM) 

-  Storekeeper  (SK) 

-  Subsistence  Specialist  (SS) 

-  Sonar  Technician  (ST) 

-  Telephone  Technician  (TT) 

-  Yeoman  (YN) 
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An  enlisted  man  (EM)  receives  an  MOS  after  either 
completing  a  specialty  school  or  by  completing  an  intensive 
on-the-job  training  progrsun.  Both  methods  are  usually 
completed  while  the  EM  is  in  the  E-3  paygrade.  An  E-3  who  has 
not  obtained  an  MOS  is  called  nonrated;  E-S's  through  E-9's 
hold  only  one  MOS  at  a  time  and  usually  keep  that  MOS 
throughout  their  careers. 

Promotions  to  the  next  highest  paygrade  are  performed  by 
MOS  and  are  dictated  by  the  needs  of  the  CG.  Promotions  to 
paygrades  E-5  through  E-9  are  vacancy  driven  and  promotions  to 
paygrades  E-3  and  E-4  occur  on  qualification  for  advancement. 
Promotion  to  paygrade  E-2  occurs  upon  completion  of  basic 
training. 

Enlisted  personnel  sign  service  obligation  contracts 
(enlistments)  that  require  the  person  to  serve  in  the  CG  for 
a  set  number  of  years.  Enlistments  usually  last  four  years. 
(There  have  been  two  year  enlistments  but  these  have  been 
terminated.)  Upon  completion  of  an  enlistment  and  upon 
approval  by  the  CG,  a  person  may  sign  a  new  contract 
(reenlist) ,  or  separate  from  the  CG.  Enlisted  attrition 
generally  occurs  because  of ; 

retirement  -  minimum  20  years  active  duty  service, 

-  non- reenlistments  -  an  EM  chooses  not  to  reenlist  or  the 
CG  chooses  not  to  reenlist  the  person, 
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-  a(3ininistrative  reasons  -  an  EM  may  depart  the  CG  prior  to 
the  end  of  his/her  enlistment  "for  the  convenience  of  the 
government" , 

-  death  or  disability, 

selection  to  an  officer  program. 

For  the  purposes  of  this  study  the  term  separation  refers  to 
a  person  leaving  the  CG  for  any  of  the  above  reasons. 

C .  PROCBDXTRB 

This  thesis  examinee,  by  paygrade,  individual  personnel 
characteristics  associated  with  attrition  and  develops  a  model 
that  calculates  monthly  attrition.  Survival  analysis 
techniques  are  used  to  study  behaviors  associated  with  these 
attrition  characteristics.  Survival  analysis  is  a  data 
analysis  technique  that  is  useful  in  examining  failure  time 
data.  Many  of  the  visible  applications  have  been  in  medical 
research  where  a  typical  "failure"  is  the  death  of  a  patient 
who  had  been  subjected  to  some  type  of  treatment.  Obviously 
it  would  not  be  practical  to  wait  for  all  the  patients  to 
"fail"  before  analyzing  the  data.  It  is  more  convenient  to 
end  the  experiment  at  some  specified  time  and  do  an  analysis 
of  all  the  data.  The  observables  that  did  not  "fail"  by  the 
end  of  the  data  collection  period  are  called  censored 


5 


observations.  They  must  be  treated  differently  since  their 
failure  times  are  yet  to  occur  and  are  hence  unknovm. 
Survival  analysis  allows  for  censored  failure  times  and 
provides  techniques  for  including  them  in  the  analysis.  The 
primary  tool  used  in  the  survival  analysis  will  b  the 
estimation  of  survivor  functions  or  curves.  These  cuzves 
provide  probabilities  that  a  person  survives  up  to  and  beyond 
a  set  time,  t.  Their  estimators  are  often  called  product 
limit  estimators  because  of  the  technique  used.  Survival 
probabilities  are  the  compliment  of  attrition  rates  and  are 
needed  to  calculate  expected  monthly  attrition  figures. 

The  personnel  data  are  cross  classified  by  paygrade  and 
each  paygrade  is  partitioned  by  certain  individual  personnel 
characteristics.  Paygrades  E-1  and  E-2  will  be  partitioned  by 
the  characteristics  of  sex,  marital  status,  and  race. 
Paygrade  E*3  will  be  partitioned  by  these  same  characteristics 
and  by  whether  the  person  had  obtained  an  MOS  or  was  nonrated. 
Paygrades  E-4  through  E-9  will  be  partitioned  by  the  same 
characteristics  as  E-ls  and  E-2s  and  by  CG  MOS.  MOSs  will  be 
grouped  according  to  categories  listed  in  Table  1  on 
page  7. 
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TABLE  1.  MOS  QROUP] 

[KGS 

Group 

Rates 

Aviation 

AD,  AE,  AM,  ASM,  AT 

Administrative 

YN,  SK 

Technical 

EM,  ET,  FT,  MST,  ST,  TT 

Skilled 

HS,  SS,  MU,  PA 

CG  Professional 

MOSS  are  grouped  for  two  reasons.  MOSs  within  the  same 
group  have  similarities  that  should  cause  them  to  exhibit 
common  attrition  behavior  and  hence  have  common  survival 
probability  values.  The  MOSs  within  each  category  all  have 
several  common  attributes;  specifically,  the  type  of  training 
needed  to  acquire  the  MOS,  the  type  of  duty  or  jobs  associated 
with  those  MOSs,  and  the  possible  civilian  employment 
opportunities  available  for  persons  with  these  MOSs.  MOSs  are 
also  grouped  to  Increase  the  personnel  Inventory  within  each 
paygrade  required  for  a  survival  analysis.  Chapter  II 
contains  an  explanation  on  these  MOS  groups  and  examines  the 
data  used  In  this  study. 

For  most  MOSs  stocks  are  large  and  empirical  survival 
probabilities  can  be  used.  Modelled  survival  probabilities 
are  needed  If  personnel  stocks  are  smaller  than  preferred. 
Small  personnel  stocks  are  subject  to  large  data  variability; 
models  are  developed  to  provide  stability.  Paygrades  E-8  and 
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E-9  could  be  considered  too  small  to  use  empirical  survival 
probabilities,  thus  regression  models  are  examined  in  Chapter 
IV  for  these  paygrades. 

A  model  that  projects  an  expected  number  (and  variance)  of 
monthly  attritions  will  be  developed  and  demonstrated  in 
Chapter  V.  This  model  is  based  on  principles  of  the  Binomial 
distribution.  Paygrades  will  be  partitioned  several  times  to 
produce  multiple  binomial  experiments  from  which  expected 
values  and  variances  can  be  calculated.  The  product  limit 
estimates  generated  by  the  survival  analysis  will  correspond 
to  the  probabilities  of  success  associated  with  a  binomial 
experiment . 

Finally  conclusions  and  recommendations  are  contained  in 
Chapter  VI.  Since  this  is  the  first  time  in  recent  years  that 
enlisted  attrition  has  been  viewed  using  operations  analysis 
tools,  most  of  the  conclusions  should  be  unique.  It  is 
expected  that  additional  studies  utilizing  more  refined  data 
will  be  needed  to  better  understand,  analyze,  and  predict  CG 
enlisted  attrition. 
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II.  DATA  OVERVIEW 


A.  POPULATION 

Data  for  this  thesis  were  provided  by  the  CG  Pay  and 
Personnel  Center  located  In  Topeka,  Kansas.  The  data  were 
provided  In  two  sets.  The  first  set  contained  Individual 
personnel  records  for  all  persons  who  separated  during  FY  91 
(1  November  1990  -  31  October  1991} .  The  second  set  contained 
Individual  personnel  records  for  all  enlisted  persons  who  were 
on  active  duty  the  last  day  of  FY  91  (31  October  1991}  .  This 
second  data  set  Is  in  a  sense  a  "snapshot"  of  the  GC  EM 
composition  on  31  October  1991.  In  both  data  sets  each 
personnel  record  contained  the  persons  sex,  MOS  (If 
applicable} ,  rank,  m&rltal  status,  ethnic  origin,  date  of 
current  enlistment ,  and  date  of  separation  (if  applicable} . 
Each  persons  length  of  service  (In  months)  can  be  calculated 
from  the  data  entries  date  of  current  enlistment  and  date  of 
separation.  It  should  be  noted  that  the  entry,  date  of 
current  enlistment,  did  not  necessarily  correspond  to  when  the 
person  first  joined  the  CO.  This  affects  the  calculation  of 
a  persons  time  In  service.  In  certain  cases  the  CG  allows  an 
EM  to  separate  from  the  CG  and  then  allows  the  person  to 
rejoin  the  CG  In  the  future  at  hie  or  her  past  paygrade.  Such 
persons'  past  service  time  Is  not  Ignored.  The  person 
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maintains  all  accumulated  service  time  for  pay  and  retirement 
purposes.  This  condition  was  not  reflected  in  the  data 
provided  for  this  study.  There  are  numerous  records  of 
senior  EMs  having  relatively  little  time  in  service  when  in 
all  likelihood  they  have  quite  a  bit  of  service  time.  This 
situation  is  most  prominent  in  the  survival  analysis  results 
of  Chapter  III.  Appendix  A  contains  procedures  and  Fortran 
progrcuns  needed  to  access  and  use  the  data. 

B.  LIMITATIOKS 

Analyzing  only  one  years  eunount  of  personnel  data  presents 
two  problems.  The  first  and  foremost  is  that  the  data  lack 
temporal  stability.  More  than  one  years  worth  of  data  is 
required  to  provide  stability.  Attrition  study  conclusions 
and  modelling  results  based  on  unstable  data  are  not  as 
reliable  as  if  they  were  based  on  time  tested  data.  Secondly, 
the  data  may  not  have  been  as  reliable  as  preferred  because 
data  inputs  could  not  be  verified.  The  FY  91  data  were  placed 
in  storage  by  the  Coast  Guard  PPC  and  were  not  accessible 
during  the  study.  Recommendations  contained  in  Chapter  VI 
address  these  shortcomings  further. 

C.  MOS  GROUPS 

As  stated  in  Chapter  I  MOS  groups  were  formed  because  of 
common  attributes  associated  within  each  group.  These 
attributes  are  listed  in  this  section. 


10 


1.  Aviation 

Individuals  with  aviation  MOSs  have  common  A- school 
characteristics  and  common  duty  stations.  (The  CG  has  three 
basic  types  of  duty  stations:  ashore,  afloat,  and  air 
stations.)  The  aviation  A-schools  (specialty  school)  waiting 
lists  are  relatively  long  compared  to  the  other  MOS  groups. 
This  lengthens  a  persons  service  time  in  paygrade  E-3,  thus 
Increasing  the  amount  of  time  needed  for  promotion  to  B-4. 
UnlDce  many  MOSs,  individuals  within  this  group  are  limited  in 
their  types  of  duty  stations.  They  spend  most  of  their  entire 
career  at  CO  air  stations. 

2.  Administrative 

These  MOSs  are  again  grouped  because  of  their  A- school 
characteristics,  choice  of  duty  locations,  and  possible 
civilian  employment  opportunities.  The  administrative  A- 
schools  have  short  waiting  lists  and  short  school  durations. 
Promotion  to  E-4  in  this  MOS  group  is  relatively  fast  compared 
to  other  groups.  Administrative  persons  also  can  serve  In 
almost  any  CG  duty  station  or  command.  Since  these  persons 
slcllls  are  mainly  clerical  In  nature  It  Is  felt  that  their 
civilian  employment  opportunities  are  limited. 

3 .  Teohnloal 

Persons  with  technical  MOSs  share  the  same  type  of 
characteristics  as  the  administrative  MOS  group,  but  with  two 
opposite  attributes.  Their  A-schools  have  long  waiting  lists. 


11 


are  long  in  duration,  and  have  low  graduation  rates  (as 
compared  to  the  other  MOS  group  A- schools) .  Like  the 
administrative  group,  persons  in  this  group  can  serve  in 
almost  any  CG  command.  However,  unlike  the  above  group, 
technically  grouped  persons  probably  have  more  civilian  job 
opportunities . 

4.  Skilled 

This  group  is  characterized  by  their  scope  of  possible 
civilian  employment  opportunities.  Persons  within  this  group 
probably  have  marketable  skills,  possibly  not  as  marketable  as 
the  technical  group  but  more  marketable  than  any  other  rate 
group . 

5.  CQ  Professional 

This  group  is  characterized  by  MOS  obtalnment  and 
employment  opportunities.  As  discussed  in  Chapter  I,  MOSs  are 
obtained  by  graduating  from  an  A- school  or  by  completing  an 
on- the  job  training  program.  This  group  is  the  only  MOS  group 
that  incorporates  on-the-job  training  programs  for  a  number  of 
their  MOSs.  MOSs  are  also  designed  for  CG  specific  functions 
(ie.  shipboard  machinists,  CG  weapons  technicians.  Shipboard 
navigators,  etc.).  These  skills  have  limited  value  for 
civilian  employment. 
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D.  SUMMARY  OP  PBRSOMNSL  STOCKS 

Table  2  provides  a  summary  of  personnel  stocks 
partitioned  by  the  rank,  sex,  race,  and  marital  status  for  all 
CO  EMs  who  were  on  active  duty  as  of  31  October  1991. 

TABLE  2.  ACTIVE  DUTY  PERSOMMEL  STOCKS  FY  91 


Total 

Male 

Female 

Cauc. 

Non- 

Cauc. 

Mar. 

Not 

Mar. 

E-l 

732 

622 

110 

595 

137 

68 

664 

E-2 

3383 

3027 

356 

2776 

607 

622 

2761 

E-3 

3651 

3358 

293 

3031 

620 

1256 

2395 

E-4 

7064 

6364 

700 

6043 

1021 

3641 

3423 

E-5 

5904 

5410 

494 

4830 

1074 

4278 

1626 

E-6 

5748 

5427 

321 

4878 

870 

4805 

943 

E-7 

2764 

2696 

68 

2436 

328 

2436 

328 

E-8 

569 

564 

5 

501 

68 

517 

52 

E-9 

289 

289 

0 

256 

33 

270 

19 

Total 

30,104 

27,757 

2347 

25,346 

4758 

17,893 

12,211 

Note  that  race  is  divided  into  two  categories:  Caucasian 
and  non- Caucasian.  There  are  relatively  few  (as  compared  to 
the  total  EM  population)  minorities  in  the  CG.  All  minorities 
were  grouped  to  Increase  the  personnel  stocks  within  the  race 
category.  As  can  be  seen  the  CG  is  a  male,  Caucasian 
dominated  service.  Table  3  on  page  14  provides  the  saune 
statistics  for  Coast  Guardsmen  who  separated  during  FY  91. 


PBESOKHSL  STOCKS  FOI 

S  SBPARATIKO  DXTKIMQ  FY  91 

Pay- 

grade 

Total 

Male 

Female 

Cauc. 

Non 

Cauc . 

Married 

Not 

Married 

E-1 

779 

637 

142 

601 

178 

83 

696 

E-2 

522 

438 

84 

429 

93 

115 

407 

E-3 

831 

731 

100 

712 

119 

276 

555 

E-4 

1248 

1094 

154 

1079 

169 

602 

646 

E-5 

438 

398 

40 

402 

36 

281 

157 

B-6 

385 

369 

16 

348 

37 

317 

68 

E-7 

315 

305 

10 

271 

44 

277 

38 

E-8 

79 

77 

2 

73 

6 

73 

6 

E-9 

46 

46 

0 

45 

1 

39 

7 

Total 

4643 

4095 

548 

3960 

683 

2063 

2850 

Approximately  the  same  proportions  of  males  to  females, 
Caucasians  to  non -Caucasians,  and  married  persons  to  persons 
not  married  separated  from  the  CG  as  did  those  who  did  not 
separate  during  FY  91. 

Table  4  on  page  15  provides  on  a  summary  of  personnel 


stocks  for  FY  91  active  duty  personnel  by  MOS  group. 


TABLE  4 .  F 

Y  91  ACTIVE 

DUTY  PBRSOMNBL  STOCKi 

3  BY  MOS  GRl 

DUP 

Paygrade 

Aviation 

Admin. 

Technical 

Skilled 

CG  Pro. 

E-3 

20 

41 

84 

71 

333 

E-4 

1142 

970 

1031 

715 

3206 

E-5 

921 

949 

781 

490 

2763 

E-6 

748 

899 

837 

527 

2737 

E-7 

262 

353 

452 

260 

1437 

E-8 

77 

88 

98 

55 

251 

E-9 

43 

40 

38 

18 

150 

Total 

3213 

3340 

3321 

2136 

10,877 

The  numbers  in  paygrade  E-3  are  small  because  not  all  E-3s 
have  an  MOS.  Overall  approximately  46%  of  all  personnel 
belong  to  the  CG  professional  MOS  group.  About  14%  belong  in 
the  aviation  group;  15%  belong  in  the  administrative  group; 
14%  belozig  in  the  technical  group;  and  9%  belong  in  the 
skilled  rate  group.  Percentages  within  each  paygrade  are 
roughly  the  same.  Table  5  on  page  16  provides  the  same 
statistics  for  personnel  who  separated  during  FY  91. 
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TABLE  5 .  F 

Y  91  SBPARI 

kTIOK  PBRSOmiBL  STOCKS 

BY  MOS  QROl 

OP 

Paygrade 

Aviation 

Admin. 

Technical 

Skilled 

CG  Pro. 

E'3 

1 

5 

12 

23 

65 

E-4 

75 

157 

180 

100 

733 

E-5 

26 

49 

91 

26 

246 

B-6 

31 

46 

51 

36 

221 

B-7 

22 

51 

23 

26 

193 

B-8 

11 

15 

14 

8 

31 

E-9 

3 

6 

3 

3 

31 

Total 

169 

328 

373 

221 

1518 

B.  DATA  PARTITIONIKQ 

The  data  was  partitioned  by  rank  to  form  9  paygrades. 
Each  paygrade  was  further  divided  by  the  personnel 
characteristics  sex«  race*  marital  status,  and  MOS  groups 
listed  in  Table  1  on  page  7. 

1.  Paygrades  B-1  and  B-2 

These  paygrades  were  divided  by  the  sex,  race,  and  the 
marital  status  of  the  individual.  Marital  status  indicated 
whether  or  not  a  person  was  currently  married. 

2.  Paygrade  B-3 

In  addition  to  the  above,  this  paygrade  was 
partitioned  by  the  obtainment  of  an  MOS.  The  person  either 
had  obtained  one  of  the  MOSs  listed  in  Chapter  I  or  had  not 
(nonrated) . 
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3.  Paygrad«i  S-4  through  B-6 

These  paygrades  were  also  partitioned  by  sex,  race, 
and  marital  status.  Additionally  they  were  partitioned  by  MOS 
groups . 

4.  Paygrades  B-7  through  B-9 

♦ 

These  paygrades  were  divided  by  race,  marital  status, 
and  MOS  group.  Sex  was  not  exeimlned  because  of  the  scarcity 
of  females  In  these  paygrades. 
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III.  SURVIVAL  ANALYSIS 


A.  INTRODUCTION 

As  Stated  in  Chapter  I,  survival  analysis  provides  a 
method  for  analyzing  censored  failure  time  data.  Failure 
times  In  this  study  refer  to  the  number  of  months  a  person 
serves  on  active  duty  prior  to  separation.  Failure  times  are 
assumed  to  be  discrete  and  are  rounded  up  to  the  next  whole 
month  figure.  For  example.  If  a  person  separates  after  31 
months,  six  days  of  service  his  or  her  failure  time  Is  32 
months.  This  chapter  reviews  the  theory  of  survival  analysis 
and  graphically  compares  survivor  functions  of  Individual 
paygrades  that  have  been  partitioned  by  single  personnel 
characteristics.  Appendix  C  contains  the  SAS  LIFBTBST 
procedure  and  Fortran  code  required  to  estimate  the  survivor 
functions . 

B.  COMPUTATION  OF  SURVIVOR  FUNCTIONS 

Let  T  be  a  nonnegative  random  variable  representing  the 
separation  time  of  an  BM  from  a  homogeneous  population.  The 
survivor  function  Is  defined  to  be  the  probability  that  T  Is 
at  least  as  great  as  some  time,  t. 
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Sit)  »Pr(Tit)  ;  0  s  t  <  «  .  (3.1) 


From  a  simple  random  sample,  this  can  be  estimated  by 


number  of  ohaervations  with  failure  time  i  2) 

total  number  of  observations 


If  T  is  a  discrete  random  variable  talcing  values  Xi  <  Xj  ...  < 
Xj,  with  associated  probability  function  f  (x^)  -  Pr(T-Xt),  then 

j 

5(t) -g  f  ,  (3.3) 

where  H  is  the  Heaviside  function,  H(x)-{0  if  x<0;  1  if  xaO} 
[Ref.  l:pp.  5-7]. 

The  presence  of  censored  data,  however,  weakens  the 
estimation  of  the  survivor  function  in  the  right  tail.  All 
that  can  be  safely  assumed  about  an  observation  censored  at 
time  t  is  that  the  unobserved  failure  time  is  greater  than  t. 
This  condition  is  relieved  by  using  an  alternate  formulation 
of  the  survivor  function  using  hazard  rates.  For  discrete 
random  variables,  the  hazard  rate,  X),  is  the  probability  that 
the  observation  will  immediately  fail  at  time  t  given  is  has 
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survived  to  time  t.  In  the  discrete  case  this  can  be 
estimated  by 


^6' 


number  of  obaezvationa  that  fail  at  t 
number  of  obaezvationa  with  failure  time  >  t 


(3.4) 


Using  hazard  rates  [Ref.  3],  the  survivor  function  may  be 
expressed  as 


5(t)-  T7  (l-X^).  (3.5) 

where  X|  is  the  hazard  rate  at  X{  for  !■!,  . . . ,  j  .  Estimating  Xj 
would  allow  an  estimation  of  the  survivor  function,  S{t). 
Chapter  2  of  Reference  1  contains  a  detailed  explanation  of 
the  computation  of  the  estimated  survivor  function. 

C.  SURVIVAL  ANALYSIS  RESULTS 

In  all  paygrades  the  marital  status  characteristic  and  sex 
characteristic  behaviors  were  identical;  married  persons 
appeared  to  have  a  higher  probability  of  survival  than  those 
who  were  not  married,  and  males  tended  to  have  higher  survival 
probabilities  than  did  females.  (See  Figures  l  and  2  on  page 
40.)  Non- Caucasians  in  paygrades  E-3  through  E-9  had  higher 
survival  probabilities  than  did  Caucasians.  Results  for  other 
characteristics  are  given  below. 
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1.  B-1 


CgLucaslans  and  non- Caucasians  have  about  the  same 
survivor  functions  until  about  two  months.  After  two  months 
Caucasians  tend  to  have  higher  survival  probabilities.  (See 
Figure  3  on  page  41.) 

2.  B-2 

Caucasians  had  higher  survival  probabilities  than  did 
non- Caucasians  up  to  about  20  months.  (See  Figure  4  on  page 
41.)  After  that  non- Caucasians  had  higher  survivor 
probabilities.  Also  note  the  significant  drop  In 
probabilities  at  24  months.  The  CG  for  a  brief  time  allowed 
two  year  enlistments.  This  drop  corresponds  to  enlisted 
persons  separating  at  the  end  of  a  2  year  enlistment. 

3.  fl-3 

Individuals  possessing  an  MOS  had  higher  survival 
probabilities  than  did  nonrates  after  22  months.  (See  Figure 
5  on  page  42 . ) 

4.  B-4  and  B-5 

The  aviation  MOS  group  In  both  paygrades  had  the 
highest  survival  probabilities  and  the  technical  group  had  the 
lowest.  The  remaining  three  groups  survival  functions  fell  In 
between  the  aviation  and  technical  survivor  functions.  (See 
Figure  6  on  page  42 . ) 


21 


5.  B-6 

All  MOS  groups  in  this  paygrade  were  very  similar  as 
shown  in  Figure  6  on  page  42.  All  had  very  high  probabilities 
of  survival  and  all  were  very  similar  in  shape.  Again  the 
aviation  group  had  the  highest  survival  probabilities. 

6.  B-7,  B-8,  and  B-9 

The  survivor  functions  in  all  three  paygrades  appeared 
fairly  flat  with  slightly  decreasing  slopes  until  about  220 
months.  (See  Figure  7  on  page  43.)  After  220  months  all 
survivor  functions  decreased  at  about  the  same  rate.  All 
survivor  functions  within  these  paygrades  are  characterized  by 
high  survival  probabilities  with  the  aviation  group  having  the 
highest  probabilities  of  survival. 
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ZV.  MODELLED  SURVIVOR  FUNCTIONS 


A.  ZHTRODUCTZOM 

RegreBBlon  modelB  are  typically  uBed  to  model  survival 
experiences  from  homogeneous  populations  [Ref.  1]  .  These 
models  allow  for  failure  times  to  be  dependent  on  observed 
explanatory  variables.  These  explanatory  variables  can  be 
either  qualitative  or  quantitative.  Two  common  regression 
models  are  the  exponential  regression  model  and  the  welbull 
regression  model.  Welbull  regression  models  will  be 
constructed  and  analyzed  for  paygrades  B-8  and  B*9.  The 
primary  reference  for  this  chapter  Is  Chapter  2  and  Chapter  4 
of  Reference  1. 

B.  WEIBULL  REQRBSSIOH  MODEL 
X.  Welbull  Distribution 

The  Welbull  distribution  Is  a  generalization  of  the 
exponential  distribution  that  allows  for  a  power  dependence, 
p,  of  the  hazard  of  time.  The  Welbull  distribution  Is  In  fact 
an  exponential  distribution  when  p»l.  The  probability  density 
function  (pdf)  o£  the  Welbull  Is 

f ( t)  fox  tiO.  (4.1) 

The  hazard  function  Is 
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x{t)  •xpat)p‘^ 


(4.2) 


and  the  survivor  function  Is 


Sit)  .  (4.3) 

A  convenient  check  to  verify  that  data  fits  the 
Weibull  distribution  is  to  plot  the  log[-log  S(t)]  versus  log 
t(d(t)  is  the  empirical  estimate  of  the  survivor  function). 
If  the  data  can  be  fitted  by  the  Weibull  distribution;  the 
plot  should  provide  a  relatively  straight  line  with  the  slope 
of  the  line  roughly  estimating  the  value  of  p. 

2.  Log-Linear  Regression  Model 

The  Weibull  distribution  can  be  generalized  to  obtain 
a  Weibull  regression  model  by  allowing  the  failure  rate  to  be 
a  function  of  a  vector  of  qualitative  personnel 
characteristics,  X.  In  this  application  X  is  a  row  vector 

consisting  of  the  binary  variables  X,vi,yoo»  . . .  X^hniwi.  X^^, 

and  X«,p^Miiotiii  which  take  on  a  value  of  1  if  the  observation 
falls  into  the  subscripted  category,  and  are  0  otherwise. 

If  the  conditional  (on  X)  hazard  is 

X(t|X)  •\pi\t)P-'e^ 
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then  the  conditional  density  of  T  Is 


f  (t|  Jf) -Xp(Xt)^-‘e*i»e-«^*''>  (4.4) 

where  |9-(/3p  . . .  ,/3a}  Is  a  vector  of  regression  pareuneters.  In 
terms  of  Y-log  T,  equation  (4,4)  Is  the  linear  model 

Y*Qi  +  XI3*+  ffW  (4.5) 

where  a- -log  X,  or-p'*,  and  W  has  the  extreme  value 

distribution.  Inequation  (4.5)  the  personnel  characteristics 
act  addltlvely  on  the  time  scale,  vice  multlpllcatlvely  In 
equation  (4.4).  [Ref.  1] 

3.  Modelled  Survivor  Function 

Equation  (4.5)  Is  known  as  an  accelerated  failure  time 
model.  It  is  characterized  by  the  scaling  of  the  failure  time 


r-e'^To 


where  Tq  Is  a  failure  time  from  a  baseline  when  the  values  of 
X  are  all  zero.  Using  this  time  scaling  the  modelled  survivor 
function  Is 
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S(t)  -e 


(4.6) . 


In  equation  (4.6),  S(t)  can  be  estimated  by  Inserting 
appropriate  values  for  X  and  t,  and  estimated  values  of  0,  X, 
and  p.  This  will  produce  a  modelled  survivor  function  which 
can  be  plotted  to  produce  curves  like  those  In  Chapter  III. 
If  this  modelled  survivor  function  Is  similar  (in  shape  and 
magnitude)  to  the  empirical  survivor  function,  then  the  model 
can  be  used.  If  It  Is  not  similar,  then  It  possibly  should 
not  be  used. 

C.  RESULTS 

The  partitioned  data  sets  were  regressed  using  the  SAS 
LIFERBG  procedure  [Ref.  4]  to  obtain  estimates  for  0,  X,  and 
p.  These  values  along  with  appropriate  values  for  X  were 
Inserted  Into  equation  (4.6)  to  estimate  modelled  survivor 
functions.  Figures  8  through  13  on  pages  43  through  46 
contain  plots  of  these  modelled  functions  against  empirical 
survivor  functions.  Generally  the  modelled  survivor  functions 
do  not  fit  the  data  as  well  as  desired. 

Nonetheless,  modelled  survivor  functions  for  these  two 
paygrades  will  be  used  pending  the  development  of  sharper 
alternatives.  The  partitioned  data  sets  are  so  small.  It  Is 
felt  that  a  model  will  produce  the  most  reliable  and  stable 
probability  estimates. 
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V.  BINOMIAL  COUNTINQ  MODEL 


A.  INTRODUCTION 

The  Binomial  Distribution  provides  a  means  for  calculating 
expected  values  and  variances  of  sums  of  random  variables  that 
have  one  of  two  outcomes:  success  or  failure.  This  study  Is 
concerned  with  estimating  the  number  of  enlisted  who  separate 
on  any  given  month  (or  conversely  the  number  who  do  not 
separate) .  This  chapter  develops  a  simple  model,  based  on 
principals  of  the  Binomial  Distribution,  that  calculates 
monthly  survival  figures  per  paygrade,  and  provides  an  example 
of  how  the  model  works.  The  primary  reference  Is  Chapter  3  of 
Reference  5. 

B.  BINOMIAL  DISTRIBUTION 

A  binomial  experiment  Is  an  experiment  that  possesses  the 
following  properties: 

-  the  experiment  consists  of  n  Identical  trials, 

-  each  trial  results  In  one  of  two  outcomes,  a  success  or 
failure, 

-  the  probability  of  success  on  a  single  trial  Is  equal  to 
p  and  remains  the  same  from  trial  to  trial, 
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trials  are  Independent,  and 


-  the  random  variable  of  Interest  Is  Y,  the  number  of 
successes  In  n  trials.  The  expected  value  of  Y  Is 
B(Y)-np,  where  n  is  the  number  of  trials  and  p  Is  their 
corresponding  probability  of  success.  [Ref.  5] 

A  person  separating  from  the  CO  meets  these  standards  and 
can  be  considered  a  Bernoulli  trial.  Each  term  trial  results 
In  one  of  two  outcomes.  A  success  Is  marked  If  the  person 
does  not  separate,  a  failure  when  the  person  leaves  the  CO. 
The  product-llmlt-estlmate  of  the  survivor  function  provides 
the  probability  of  a  person  surviving  to  and  beyond  a  time,  t; 
P(Tiit)  for  t-l,...oo.  This  probability  estimate  remains  the 
8aii\e  for  each  trial  (le.  for  all  persons  with  length  of 
services  equal  to  time  t)  .  The  decision  by  a  person  to 
separate  (or  the  COs  decision  to  separate  a  person)  Is  not 
affected  by  another  persons  choice  and  thus  independent  of  any 
other  persons  decision.  And  finally,  Y  equals  the  number  of 
persons  not  separating  within  a  paygrade  who  have  survived 
beyond  a  time  t.  Using  these  definitions  a  basic  model  can  be 
constructed  that  generates  an  expected  value  and  variance  for 
the  number  of  persons  that  do  not  separate  on  a  given  month. 
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C.  COUMTINQ  MODEL 

Paygrades  are  partitioned  in  two  steps  to  form  data  cells 
that  consist  of  Independent  and  Identical  trials  required  of 
the  Binomial  Distribution.  The  first  step  separates  each 
paygrade  Into  the  single  personnel  characteristic  examined  In 
the  survival  analysis.  These  will  be  called  sets.  The  second 
step  further  partitions  the  paygrade  by  the  eunount  of  time 
each  person  has  served,  or  Time  In  Service  (TIS)  cells. 
Figure  14  provides  a  tree  diagram  exeunple  of  this  partitioning 
scheme. 


Figure  14.  Counting  Model  Partitioning  Soheme 


T  is  the  maximum  TIS  any  one  person  has  within  each  set . 
Paygrades  E-l  and  E-2  are  partitioned  by  race  to  form  two  data 
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sets  each;  paygrade  E>3  Is  partitioned  by  MOS  obtainment;  and 
paygrades  E-4  through  E-9  are  partitioned  by  MOS  group.  Each 
cell  has  a  corresponding  probability  of  success,  pj,  (either 
estimated  from  the  actual  data  or  generated  by  a  model)  which 
wa&  defined  in  Section  4.D.1  to  be  the  probability  of  survival 
being  greater  than  or  equal  to  t.  Given  nj,  the  expected 
number  of  persons  who  fall  into  cell  j ,  the  expected  number  of 
successes  (persons  that  do  not  separate)  in  each  cell  is 

E(Yj)'‘a^Pj  ;  (4.7) 

The  cell  variance  is 

VariYj)  (4.8) 


where  gj-l>pj.  The  expected  number  of  successes  within  each 
set  is  the  sum  of  the  cell  expected  values 


r 

E(Y,^)  •J^EiYj) 


(4.9) 
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and  by  svunmlng  cell  variances  the  set  variance  Is 


T 


Var(YJ  VariYj) 
;■» 


(4.10) 


The  expected  nxunber  of  successes  per  paygrade  Is  calculated  by 
summing  set  expected  values.  Paygrade  variance  Is  the  sum  of 
all  set  variances. 

Appendices  C  and  D  contain  procedures  for  generating 
product  limit  estimates  which  correspond  to  the  pyS. 

D.  EXAMPLE 

Data  from  paygrade  E-1  is  used  to  illustrate  this  counting 
model.  Paygrade  B-1  contains  2  paygrade  sets.  Set)  contains 
records  for  Caucasians  and  Set]  contains  records  for  non* 
Caucasians.  Table  6  and  Table  7  on  page  32  contain  frequency 
counts  of  the  number  of  persons  falling  within  each  cell  per 
set. 
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Tables  8  and  Table  9  on  page  33  contain  the  estimated  pj  values 
for  each  data  sets. 
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TABLE  9.  PRODUCT  LIMIT  BSTIMATBS  PAYQRADB  B>1  MOK-CAUCASIAN 


Substituting  hj  the  values  from  tables  6  and  7  for  the  nj 
values  and  f>j  the  estimated  probability  values  from  Tables  8 
and  9  into  equation  (4.9),  the  estimated  expected  value  of 
each  set  is 
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-353.475 


iHYi) 


a 

-524^>“7S.905 


where  eet|  corresponds  to  Caucasians  and  set}  corresponds  to 
non- Caucasians .  Using  equation  (4.10)  the  ^ar(Y,)  -  129.5  and 
Var(Yj)  -  30.97.  The  overall  paygrade  estimated  expected 
value  and  variance  is  the  sum  of  the  set  values 

=430.38 

(-1 


For  working  purposes  these  values  are  rounded  down. 
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VI.  CONCLUSIONS  AND  RBCOMKBNDATIONS 


A.  CONCLUSIONS 

There  are  approximately  30,100  enlisted  persons  In  the  CG 
at  any  given  time.  A  method  Is  needed  to  predict  monthly 
attrition  rates  to  maintain  the  CG  enlisted  hierarchical 
structure.  Survival  analysis  was  used  to  examine  the 
characteristics  of  sex,  marital  status,  race,  and  MOS  group 
behavior  of  CG  enlisted  personnel.  The  results  were  fairly 
interesting.  Males  and  married  persons  tended  to  have  higher 
probabilities  of  survival  than  did  females  and  individuals  who 
were  not  married.  Non- Caucasians  in  most  of  the  paygrades  had 
higher  survival  probabilities  than  did  Caucasians.  MOS  groups 
demonstrated  different  characteristic  behaviors  among 
different  paygrades.  MOS  groups  In  paygrades  E-4  and  E-5 
behaved  about  the  same  up  to  48  months.  After  48  months  all 
groups  showed  a  rapid  decline  In  survivor  probabilities  with 
the  aviation  group  having  the  highest  probabilities  and  the 
technical  group  having  the  lowest.  This  drop  Is  Indicative  of 
the  end  of  a  four  year  enlistment  when  Individuals  can  usually 
separate  from  the  CG.  MOS  groups  In  paygrades  E-6  through  E-9 
all  had  very  high  probabilities  of  survival  and  demonstrated 
a  decline  at  certain  months.  E-6  showed  this  decline  at  about 
200  months  and  E-7  through  E-9  demonstrated  this  at  about  240 
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months.  The  drop  in  probabilities  at  240  months  is  easily 
explained.  Many  Coast  Guardsmen  retire  at  20  years  and 
receive  a  substantial  pension  and  many  important  benefits  (ie. 
medical  coverage  for  the  person  and  his  or  her  dependents) . 
The  drop  at  200  months  In  paygrade  B-7  is  an  anomoly.  An  E-7 
at  this  point  in  his/her  career  has  almost  17  years  of  service 
time  and  has  only  about  three  more  years  to  serve  before 
retirement.  Separating  at  17  years  seems  to  be  financially 
inappropriate  because  the  person  will  not  receive  a  pension 
nor  benefits  if  they  separate. 

Modelled  survivor  functions  were  developed  for  paygrades 
E-S  and  E-9  because  these  paygrades  possessed  small  population 
sizes.  Despite  the  fact  that  the  models  produced  generally 
poor  fits,  it  was  felt  that  these  survivor  functions  would  be 
more  reliable  than  survivor  functions  generated  from  very 
small  data  sizes.  Therefore  models  were  used  to  calculate 
attrition  in  paygrades  E-8  and  E-9. 

The  binomial  count ixig  model  provides  an  effective  tool  for 
calculating  monthly  attrition  figures.  Expected  values  and 
variances  can  be  calculated  by  partitioning  each  paygrade  into 
several  binomial  experiments.  This  model  is  easy  to  use  and 
can  be  expanded  in  future  studies. 

B.  RECOMMENDATIONS 

As  a  first  study  the  goals  of  this  thesis  were  met. 
However  further  studies  are  required  to  better  analyze  CG 
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enlisted  attrition.  The  counting  model  could  be  adjusted  or 
refined  based  on  findings  from  the  recommendations  listed  in 
this  section. 

1 .  Recommendation  1 

The  counting  model  needs  to  be  verified.  Personnel 
data  from  early  months  of  FY  92  need  to  be  processed  through 
the  model  to  verify  that  it  provides  better  results  than  the 
current  calculation  method. 

2 .  Recommendation  2 

Additional  historical  data  needs  to  be  analyzed  to 
provide  temporal  stability  to  the  study.  The  counting  model 
produced  in  this  study  assumes  that  the  characteristic 
behaviors  of  individuals  identified  in  FY  91  will  repeat 
themselves  in  FY  92 . 

3 .  Recommendation  3 

A  more  manageable  data  base  needs  to  be  Identified 
that  allows  for  accessible  man-machine  interface.  An 
accessible  data  base  allows  the  analyst  to  verify,  manipulate, 
and  easily  examine  data  entries  as  an  on-going  task.  This 
provides  data  reliability.  The  CG  PPC  is  not  responsible  for 
maintaining  this  type  of  data  base,  thus  PWP  should  examine 
this  issue. 

4 .  Recommendation  4 

Individual  MOSs  within  MOS  groups  need  to  be  exeunined 
to  ensure  they  demonstrate  similar  attrition  behaviors. 
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5.  Raaommandatlon  5 

Combined  personnel  characteristics  need  to  be  analyzed 
for  survival  behaviors  and  seasonal  attrition  trends  need  to 
be  exeunined. 

6 .  Reaonmendation  € 

The  data  base  needs  to  be  extended  In  time.  At 
present  there  Is  no  justification  for  assuming  time 
homogeneity  of  rates.  Such  an  extension  would  also  allow  the 
study  of  some  important  exogenous  variables,  such  as 
unemployment  rate  and  Inflation. 

7 .  Reooamendation  7 

The  educational  levels  of  Coast  Guardsmen  should  be 
analyzed  for  their  effect  on  attrition  behaviors.  High  school 
graduates,  alternate  high  school  degree  holders,  and  non -high 
school  graduates  need  to  be  considered  to  determine  if  there 
are  differences  in  attrition  behaviors  amongst  these  groups. 
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APPENDIX  A.  PROCBDUABS  FOR  ACCBSSINQ  DATA 

1.  Data  is  stored  on  two  magnetic  disks.  To  access  disks: 

a.  login  to  mainframe  computer 

b.  type  "getmvs  mv8215  (separation  data  tape) "  or 
"getmvs  inv82l7"  (active  duty  data  tape) 

2.  To  run  a  program  that  has  a  huge  output  file,  ie.  a  file 
that  will  not  fit  on  an  A-disk: 

a.  type  "tdisk  6"  to  get  a  temporary  disk  with  6 
cylinders 

b.  switch  A'disk  with  temporary  disk  (  say  temp  disk  was 
assigned  "D  130" 

1.  type  "access  130  A" 

2.  type  "access  191  D” 

c.  copy  program  from  the  now  D-disk  to  the  newly  created 
A-disk 

d.  run  program 

e.  switch  disks  back  when  finished  by  reversing  procedure 
in  2.b 
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APPENDIX  B.  FiaURBS 
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Figure  3.  Survivor  Function  Faygrade  B-1  by  Raoa 


(•Avaiuoe  E-a  by  race 


Flgvra  4.  Survivor  Function  Paygrada  S-2  by  Race 
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Vlguro  8.  Modal  Pit  Faygrad*  B*8 
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Plgur*  13.  Model  Pit  Paygrad*  B‘9 
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APPENDIX  C 


PROCEDURES  FOR  OBNERATINfl  SURVIVAL  FUNCTIONS 


rikCi  LIPCtItT  CSKI  *1 


ly  DMI'  /•  THK  l«l  eOH  TAKIt  A  BATA  Ut  AND  CtNMATII 
lUAVtVOA  AUNCItONt  ADA  THI  DATA.  THIt  COBt  HAt  MID  ADA  AAVDAADI  l-l 
AND  eilATII  IMVIVAC  AUNCTIOM  AOA  CAUCAIIANI  and  NOW-CAUCAItANI,l/| 

BATA  AOil 

INAUT  AATI  AAM(  tIA  ABC  HU  lACI  HONTHI  CIH  /aTNII  LIHI  DfAINCI  THI 
yAAIAMI  HlTNlN  THI  DATA  UT.A/I 

BAOUAilAACI<ll  /•THII  tIAUATCt  AkL  CAUCAtlANA  AAOH  NaH>eAUCAtlANt . 

I  InDICATII  DAJCAttAN.  THIt  lAUlTl  THl  DATA  UT  INTO  1  IITl  1  CDNTIANI 
CAUCAIIANI,  THI  OTHIA  CONIAlNI  AAL  DTHCA  NINOAITITII.a/i 

CAABI  /ATILAI  lAI  THAT  THC  BATA  Uf  DOLLOUt.  DATA  ilT  NUIT  II  ALACIB 
IMHEDIATIUV  AATtI  THil  UNl  AATIA  THI  UMICOCBN,  AND  HUlT  HAVI  A 
IIHUOLON  ALACIO  IHHIBIATILV  AATCB  THI  LAIT  Bill,*/  I 


Blue  LIBITIIT  BLOTiIII  MAKTI«<AI  OUTIIAVUOCIIT  />TH1I  II  THI  MAIN  CODI 
II  CKAltl  A  tUlylVAU  BUNCTION  BkOT  Itl.  AND  klMITl  TH|  N-ANII  DA  iHl 

Blot  to  «i  monthi.  aiso  tiou  thi  bioduct  limit  iitinatii  in  a  sas  bui 

CALlID  BDCCIT.A/  I 

riMI  MOMHIiClNlll  /IDIBINII  TW  T|Mt  VMlAiLt  IMONTHtI,  AND  lITt  Thi 
CINICIINB  VAIIaIlI  HITHIN  TM  DATA  UT  COUAL  TO  I  •/  I 

ITIATA  BIOUB  /aBIABOIMt  THI  BIOCCDURI  ON  ALL  DtBIHID  OlOUBt,  IN  THIS  CAU 
TnIII  Alt  ONLV  I  OIBUBIA/  I 

Bloc  BUNT  OATA>BOillT  /iBIINTt  Tl€  DATA  UT  BOUIT  INTO  A  CMI  Bll.li/I 


BUI,  LIBIIIIT  CODI  At 


TITLI'BAVOIADI  l-T  BY  lATI  BIOM  '  /UAMI  At  BIIVIOM  CODI  INCIBT 

IT  OlViDtl  THI  BATOIADl  INTO  I  OlOUBII  I  BOA  l‘CH  MOl  OlOUB  .  BIUUB  I 

ID  THI  AYlAIION  BIOBLt,  liAOMIN  ■  liTCCMNICAL >  AiWKLID'l'CO  BIO  •/  1 

'  , 

data  BOTi 

IHBUT  lATI  bank  1(1  ACI  HAI  lACI  MONTHI  CIMl 
IB  l,AriiUOIIAt|,17tllAUitl»IIAT|t||OIAATI«ll«  TMIMOBill 
ILCt  IB  llt(•l(0|BA1l•A|l  THIN  SBill 

tl.ee  IB  IATIi)IOIBAUi|(OIIATI»lABI»ATI«lUIIATliJ40IIATliSIO  THIN  0Bi|l 
iLtl  IB  IAU*I0(IIA1I.|40IIATI>BU  TWN  OBiAl 
(I.U  OB  >11 
CAlDtl 

I 

BIOC  LIBITIIT  BLOTHII  HAKTIMlilBO  OUTUIV'BOUITi 
TIHl  MOhTHB>C(NII ll 
ITIATA  BBi 

BNUC  Bllhl  DATAiBOTIITl 
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INTCtlt  I.M.N.r 


wdfit; 

THII  FII03FAN  OtHCaATCI  MaC[t.l.tO  raODUC*  LIMIT  ||T|HAtl|  FM 
FArOFADt*  HMICH  HAVf  IttM  FA«T|TI«MCD  IV  I  mt  OFOUFtI  AVU1ION. 
ACATlMieTKAlIVt.  riCHNICALi  SKlLLtC'  CO  F*0*IS>IWIAL.  IT  Kli  M 
IMFUf  UtOaiUieM  COtFFICIINTt  IKTAl  TMISUOH  ilTAtl).  tl«AFt 
FAAAMfnat  IF  I  AND  Fil.  tCALI  FAFAMITfFI  ILAMIDA)  AMS  LAMtM)< 

AND  A  TIMC  FAHAMdll  IH.  tMli  FFOUAAM  It  ALSO  OtttONIt  TO  HOOX 

ON  data  that  fit  a  I  Fiiei  miiial.  vaaiailis  ktai  ttav  kta4< 

FI.  AND  LAHIDAl  aFI  VALWI  FOO  TMI  FUST  MliULL  FIICI  IF  It  TNI 
TIHC  yALlX  OF  IMAI  THI  FIFIT  iFCAF  OCCI^tl.  yAAIAILII  NTAt  TItU 
HTAO.  FI.  AND  LAHIDAl  IAN»  Ml  COFFItFONS  TO  ?»«  tICOHO  FIta. 
VAOUILtS  tITAF  THOU  IITAII.  FI.  LAHIDAl  lANO  Ml  COMtlFWO  TO  TT« 
last  NIIIVA.L  Mliel.  VAAIAILI  UTAH  ICUALt  IIFO  AND  COMtlFONn 
TO  THI  CO  FFOFtStlONAL  eOOFFtCItNT.  TMI  VAFIAILt  tMAT  It  TIC 
HODCLlID  FAODUCT  LIMIT  ISTIMATIl  IMATI  FOA  TIC  AVIATION  Hit 
OAOLF.  IMATI  FOA  ADMIN  OAOUF,  IMATI  FOA  TICTWICAL  OFOUF.  ITUTA  FOR 
THI  INILLID  OAOUF.  AND  IMATI  FOA  TIC  CO  FAO  OAOLF.  SHAT  VALWI  AR 
CALCULATCO  FOA  TIMt  FCAIODIT.  (Til .... .  NT.  UCAI  N  lA  TIC  LAIT 
TllC  VALUl  MtlDID  FOA  TMI  FAYOAAOt, 

THI  IITIMATtI  AAI  NAITTIN  ON  THI  CHI  FILIII  FTIIFIOI.  FTOIFIOli 
AND  FTTDFtOI , 


\  FCAL  LAHFClM'MCTAl.KTAt.rtTAl.lCTAA.FtTAt.lCTAIl 

FCAL  LAMHA.VFr.lt'AI. IITAMCTAl.ItTAI. IITAII. FtTAlMITAII  ■ 
AfAL  IMATI. CMATi.lHATMMATA. IMATI. LAMIOAtiFI 

,  C  FOLLONINO  ItTAFLIINCI  ALL  VMIAlLt  VALWI. 

FAAANITIAINi|O.LAM|DA|a.|eiT  .F|it.AIII 
FAAAMtTf A  lit TAlie.I.IITAli-, lilt. ltrA|i.,»t||.|tTA4i>  ,11111 
FAAAN|T|AIN<ICe  .LAMIDAltO.IDOOl  .Flll.ilSI 
'  FAAAICTIAIIITA|i-.«IA|.|tTAti.llf.lirAT>-,»IIT,|tTAO»4.TIT«l 

FAlAT>tTtAIFiA]«>LAM|gAl<.ie;i  .FUD.OIH  ) 
FAAAMtTIAIIITA|i.e;iA.KTAIDF.II>MltA||<.|l«tllTA|l>-,lllll 
FAAAMttlAIIITAIWe.l  I 

e  CALCAATII  IMATS  FDA  FIAIT  WIIUtL  F|ir|, 

DO  II  IM.N 

IHATIlCaFI-TIlAMIDAlAliINFI-DITAIlliiFII) 
IMAr:i|aF|>tlLAHIJAl>l>tllAI-IITA3)|llFIII 
'■  tMATIiI«F(-IILAMIDA|l|iIIIFI-D|TAII|MFIII 

INATAi|l(F|.|ILAMiDA|>|ltKAI>ltTAI))>>FIII 
tHATtit>FI>IILAMIDA|i|>IXrl-DtTAIII)l>FIII 
NA nl III . II 11  . IHAT I , IHATI , IMATI. IHATA .IMATI 
IF  COMTINUt 

C  CALCULATII  SMATS  FOA  CICDMD  NItAULL  FIICI. 


DO  it  l>M>l .N 

|HATI>(XF|.|ILAM|DA:i|l|)(FI-DITA|)|liFtl| 
tMATi>CAF|.ilLAHtOA;<l>t«Ft-ICTA|)|TtFil) 
tMATliCXFI-l|I.AM|CA;M»tAFI-ltTAn»>F;)) 
IHATA<IXFj.|ILIMIDA;>l>tAFI>ltTA|))l>F|)) 
tMAT|>IXF|-IILAM|eA:>|ilAFI>ltTA|!l)>>Fi)l 
MAITIIIA.TI tl.CMAII  .IHATi.lMATI.IMATA.IHATI 
it  eONTIMUl 

C  CALCULATII  IHATI  FOA  LAIT  FIICI. 


DO  It  liN.I.F 

SHAT  I  >(XF  |.  I  ILAHIOAnlUAFI-IITAt  ll■■FI  I) 
SHATt.|XF|.|ILAMIDA|«l>II(FI-IITAI«ll>>l>ll) 
IMAT|itX'(-(ILAM|DA|<|i|IIFI-IITAIIIIliF|ll 
SHiTA.[XFI.|ILAMICA|x|itXFI-D|TA||i;i>.'|)) 
SHATS>[XF|. I>LAMICA|>|iCXFI-ICTAIIITx>F))I 
MFITtITO.TI II iSHATI.  IHAT], SMATI.SHATA, IHATI 
It  eoNTINUI 

II  rOAMArilX.ll.iX.FA.A.IA.FI.A.IA.Ft.A.IX.FI.I.lX.FA.AI 
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APPENDIX  D.  PROCEDXTRSS  FOR  CONDUCTING  REGRESSION  ANALYSIS 


fiooRAH  mor«ir 


THU  CODt  ItTADtTII  A  AAVOKADI  DAf*  >IT  IV  HOI  OAOl»  ANO  CIIAfll 
BATA  tltt  NCIBID  TO  CONOUCT  THI  lAI  IIAIAEB  AVKIDUAI.  AOAIN  THU 
CODC  It  BtSieutD  AOA  A  I  IIICI  KIIIA.L.  THt  IHAU1  TUI  Ull 
CONTAINS  THt  AAVOAtDf  tiATA  tfT,  1  OUTTUT  OATA  tltt  AAt  CAtATIB  TO 
CONTQAH  TO  THt  I  AltCK  HtllUI.1  .  THI  VAAIAlLC  N  It  THI  NUHKA  OH 
HILII  IN  THT  AAVOAADI. 

IHTTOIA  I.NiAATt.HONTHSiCCH 
AAAAMITtAIHi«AII 
OTtNIi.f  aO'AOi  DATA  A') 

DO  10  lihN 

AlAOUitSI  AATIiMONTHtiCIN 
THU  It  IHt  AIAST  met  Of  THt  HtllAl. 

If  IHONTHt  .Ul,  It  I  THCN 

IflAATt  .10.  t:0  ,0A.  AATC  ,10.  I«0  .OA.  AATt  .lO.tTO  .0*. 

I  AATt  .ITI.  lit  .OA.  AATI  ,10.  ItOITHCN 

HAITtlU'.tn  HONTHSiCIN 
ILtriflAATt  .(0.  tit  .OA.  AATt  .10.  Attl  THCN 
HAITtIU.tAI  HOHTIIt.etH 

CLSCIflAATI  .10.  ITO  .OA.  AATt  .CO.  TAt  .OA.  AATt  .10.  Ill 
I  .OA,  AATI  ,tO.  Tie  .OA,  AATt  .CO.  lAO  .OA.  AATt  .tO. 

I  :ttl  THCN 

HIITICA.TII  HDNTHt.CCH 

CLtCIAIAAII  .CO.  ITO  .OA.  AATI  .10.  lot  .CO.  AATI  ,10.  lAOITHIN 
"AltlCA.tll  HOHTHI.OIH 

CLil 

HAITIIil.tll  NOHTHt.CtN 
INDIA 

THU  It  Tw  tceoHD  fiiei  of  thc  hiiiull. 

ILtCIf  INCHTHS  .OT.  10  .AND.  MONTHS  ,1.1,  iOOITHtH 

IflAATI  .10.  i;i  ,0A.  AATI  .10.  lit  .OA.  AATt  .10.111  .OA. 

I  AATI  .10.  ill  .OA.  AATI  .10.  IIOITHIN 

MAlTIItT.tll  HOHTHt.CIH 
lUStlflAATI  .10.  lAO  .OA.  AATI  .10.  <901  THtH 
HAITICT.tAI  MONTHS. CCH 

ILSIlflAATI  ,10.  Ill  .OA.  AATI  .CO.  JAO  .OA,  AATI  .10.  Ill 
I  .OA,  AATI  .10.  TOO  .OA.  AATI  .10.  UP  .0*.  AATI  .10. 

1  ItOl  THCN 

UIITCi:9,ll)  MONTHS. CCH 

IttllKAATI  .10.  lie  .OA.  AATt  .tO.  tOO  .OA.  AATI  .10.  KOITHIN 
HAITIISI.SII  MONTHS, CCN 
iLse 

Hinilll.TII  MONTHS, CtN 
tNDIf 


e  THU  U  THI  LAST  fllCI  Of  THI  N|lltA.t.. 

rti 

IflAATI  ,10,  Ite  .OA,  AATI  .to.  lit  .OA.  AATt  .10.110  .OA. 

I  AATI  ,10,  III  .OA.  AATI  ,10.  IIOITHtN 

NAITIlUitll  MONTHS. CIN 
iLirtflAATI  .to.  Ite  .OA,  AATI  .10.  Alei  THIN 
NAITICl.UI  MONTHI.CIN 

ILSIlflAATI  ,10.  lie  .OA.  AATt  ,10.  ttO  .OA.  AATI  .10.  110 
I  .rt.  AATI  ,10.  Tte  .OA.  AATI  .10.  lie  .oa.  aati  .10. 

1  lie)  thin 

NAITIlII.etl  MONTHS. CIN 

ILtCtflAATI  .10.  ITS  .OA.  AATI  .10.  101  -OA.  AATI  .10.  KOITHIN 
NAITFtrA.eil  HONTHSiCIN 

ILII 

NAITIIIS.SM  HONTHSiCIN 
INDIA 
INDIA 
IS  CONTINUI 

C  THU  IS  THt  fOAHAT  FOA  THI  AIAD^INfUT  STATIHINTS, 


09  fOAMATIIAiliillAilliU'lll 

THItI  ACAMAT  STATIHINTS  CAIATI  AtOAIISieW  DATA  SITS.  HONTHI  IS  Til 
ptflHDINT  VAAIAILl  AND  It  lOCAlIB  IN  THt  flAIT  CCtOMN.  THI 
AAICIDINO  columns  CONTAIN  A  0  OA  I  lA  THI  ‘J’® 

I••rf  THt  HOI  OAOOfl.  (THt  STH  HOS  OAOUf  11  ASIVMtD  ANO  DBlt  NOT  SHO 
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